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Abstract 

The term "Active Plasma Resonance Spectroscopy" refers to a class of diagnostic methods which 
employ the ability of plasmas to resonate on or near the plasma frequency. The basic idea dates 
back to the early days of discharge physics: An signal in the GHz range is coupled to the plasma 
via an electrical probe; the spectral response is recorded, and then evaluated with a mathematical 
model to obtain information on the electron density and other plasma parameters. In recent years, 
the concept has found renewed interest as a basis of industry compatible plasma diagnostics. This 
paper analyzes the diagnostics technique in terms of a general description based on functional 
analytic (or Hilbert Space) methods which hold for arbitrary probe geometries. It is shown that 
the response function of the plasma-probe system can be expressed as a matrix element of the 
resolvent of an appropriately defined dynamical operator. A specialization of the formalism for a 
symmetric probe desing is given, as well as an interpreation in terms of a lumped circuit model 
consisting of series resonators. We present ideas for an optimized probe design based on geometric 
and electrical symmetry. 



I. INTRODUCTION 



Of the many diagnostic techniques available or proposed for low temperature plasmas, 
only a few are compatible with industrial requirements. A diagnostic tool which is useful for 
supervision and/or control of technological plasma processes must be i) robust and stable, 
ii) insensitive against perturbation by the process, hi) itself not perturbing to the process, 
iv) clearly and easily interpretable without the need of calibration, v) compliant with the 
requirements of process integration, and, last not least, vi) economical in terms of investment, 
footprint, and maintenance. 

In this paper we focus on a one very promising approach to industry compatible plasma 
diagnostics, the so-called plasma resonance spectroscopy. This approach attempts to ex- 
ploit the natural ability of plasmas to resonate on or near the electron plasma frequency 
u pc = e 2 n c /e m c for diagnostics purposes. The idea to use these resonances as a measure 
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and has found 



for plasma parameters dates back to the early days of discharge physics 
many different realizations since then 2Hl5||. We follow the classification of 16[ to structure 
the many different realizations: Active plasma resonance spectroscopy couples a suitable 
RF signal into the plasma and records the response 2Hl3j|. Passive plasma resonance spec- 
troscopy merely observes the pre-existing excitation of the plasma by other sources (typically, 
the RF power input) Q, [l^. Electromagnetic concepts are based on the interaction of the 
plasma with the full electromagnetic field, i.e., electromagnetic waves. They operate above 
the plasma frequency and observe the plasma-induced shift of transmission line or cavity 
resonances which are already present in vacuum {2-5]. Electrostatic concepts, on the other 
hand, result from the interaction of the plasma with the electric field alone; they utilize res- 
onances below Up C which are not present in vacuum (standing surfaces waves) and deduce 



. Finally, one-port 



the plasma parameters from the absolute value of the frequency 

concepts utilize probes with only one electrical input/output port |4|, |5|, I7l-ll5|. while two- 
port concepts use a probe configuration with two ports or even two separate probes {2 6]. 
(Multi-port concepts would also be conceivable.) 



Numerous works deal with simulation and modeling of specialized aspects " 
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such diag- 
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22|. In 



nostic devices to achieve a better understanding of the resonance behavior [6|, |12|, 
contrast, the purpose of this manuscript is to study the whole class of active, electrostatic 
methods in terms of a general description based on functional analytic methods in order to 
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gain a deeper understanding, as no particular assumptions on the probe geometries or the 
homogeneity of the plasma are made. 



II. A MODEL OF AN ELECTROSTATIC PROBE IN ARBITRARY GEOMETRY 

We represent the plasma chamber (see figure [TJ as a simply connected, spatially bounded 
domain V, most of which is plasma (a simply connected subdomain P). Other subdomains 
of V are the plasma boundary sheath S, which shields the plasma from all material objects, 
and possibly dielectric domains D. The boundary dV of the domain V is either grounded (G) 
or ideally insulating (/) with vanishing conductivity and permittivity. 

Into this idealized plasma chamber, a diagnostic device is immersed. To render our con- 
siderations general, we assume that it is a probe of arbitrary shape, containing an arbitrary 
number of powered electrodes. The electrodes surfaces Ek, k = 1 . . . N are also seen as part 
of the domain boundary dV, of course they are insulated from each other and from ground. 
A possible dielectric shielding of the probe is represented as a part of the subdomain D 
within the plasma chamber V. Fig. 1 illustrates the assumed geometry. 

Within the subdomain P, the plasma is described by the cold plasma model, a set of 
two coupled linear equations for the RF charge density p and the RF current density j. 
The first expresses the conservation of charge, the second describes the acceleration of the 
electron fluid by the RF electric field E = — V$ and the momentum loss due to electron- 
neutral collisions measured by its rate z/, 



The electric potential $ is not only defined in the plasma domain P, but also in the sheath S 
and in the dielectric D. RF charges are present in the plasma P and - as surface charges - 
on the plasma boundary dP. The sheath, however, is electron-free and hence RF charge free. 
Consequently, Poisson's equation gives the RF potential as 




(1) 



(2) 




p in the plasma P and on dP, 



in dielectric domains D and the sheath S. 



(3) 
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The respective permittivity of the domains is given by 

{1 in the plasma P and the sheath S 1 , 
e r (r) in dielectric domains D. 

The electrodes Ey. are driven by RF voltages Vy, with k — 1 . . . N. On the grounded parts 
of the wall or the probes the potential vanishes. It is convenient to treat these sections as a 
further electrode Eq with a vanishing applied voltage Vq = 0. Taking also into account the 
insulating surfaces /, the boundary conditions for the potential are of mixed type: 

Electrodes E k : $| Bfe = V k , k = 0...N 

(5) 

Insulators /: n-V$|/ = 0. 

The currents carried by the electrodes are the measured signals of the diagnostic. 
As the electrodes are shielded from the plasma, at least by the sheath, the currents are 
carried by displacement alone. Counting them positive when they flow from the electrode 
to the plasma (in the direction of the negative surface normal), one obtains 

h = / eo^r^— • a r. (6) 
J Ek at 

Of course, there is a standard approach these equations: Applying the Fourier transform 
to the equation of motion leads to the RF version of Ohms law, an expression of the total 
current in terms of the field. The model reduces to an elliptic equation in the domain V, 

V-(e e (off) V$) =0, (7) 

with the effective permittivity e( cff ) given by 

,2 



e (eff) 



1 ? — pe . \ in the plasma P, 

00(00 — %V) 



1 in the sheath S, (8) 

e r (f) in dielectric domains D. 

Together with the Fourier transformed boundary conditions (jSJ), the elliptic equation can 
easily be solved numerically. Its linearity allows to write the solution $(r) as a superposition 
of N fundamental solutions, 

N 

$(f) = ^^(f). (9) 
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each of which obeys equation (JTj) and the boundary conditions (Ski is Kronecker's delta) 

Electrode E k : ^\ Ek = S kh l = 0...N, 
Insulators: n-V$i|ar = 0. 

With the help of the fundamental solutions a response matrix can be defined as 
the flux of the Fourier transformed displacement current through electrode Ek, 

Y M = [ iue e^V^-d 2 f. (11) 

JE k 

The Fourier transformed RF current Ik as the probe signal can then be calculated as 

N 

h = J2 YkiVi >- ( 12 ) 
1=1 

This approach has been successfully employed for numerical investigations of the plasma 
absorption probe and the multipole resonance probe in [l^] and It relies, however, on 
numerical calculations and gives only a limited insight into physics of the probe response: 
Although the mathematical structure of © to (fT2l) is quite transparent, it is not possible to 
separate the plasma contributions from the influences of the geometry from the numerical 
solution. 



III. CURRENT AND ENERGY BALANCE EQUATIONS 

With the previous section the model is defined: We now focus on the balance equations 
of current and energy and prove the physical correctness of the model. We are finally able 
to set up a general and abstract model. 

In order to obtain a deeper understanding of the probe response, we continue here with a 
general formalism. We split the electric potential into two parts, one of which is the vacuum 
potential which can be written as a linear superposition of characteristic functions, 

N 

$(f) = <f)(f) + $( vac )(f) = 0(r) + V i^i(r) ■ (13) 

1=0 

The characteristic functions with I = . . . N depend only on the geometry of the chamber 
and the geometry of the probe. They are independent of the plasma and its spatial distribu- 
tion. To calculate the the related Laplace equation has to be solved, i.e. V-(eoe r V\l/;) = 0. 
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The boundary conditions are given by 

Electrode E k : *,|^ = 6 kh l = 0...N, 
Insulators: fi-V^lar = 0. 

The first part of the electric potential defined in ( Tl3|) . referred as the inner potential 0, is 
given by Poisson's equation in the plasma and the sheath, or dielectric zones, and satisfies 
homogeneous boundary conditions. One obtains the boundary value problem 

p in the plasma P and on dP, 
-V-(e o e r V0) = { (15) 
in dielectric domains D and the sheath S, 

Electrode £ fc : (j)\ Ek = 0, k = 0...N, 

( 16 ) 

Insulator: n-V(p\di = 0. 

With the knowledge of the two distinct parts of the potential, we are able to reformulate 
dynamical equations for the charge density and the current density. From ([T|) and (j2J) we 
ultimately arrive at 

|-V7. (17, 
d* N 

-i = -e O wJeV0 - i/J - e c^ o ^ ^V^ ? . (18) 

1=0 

Similarly to the potential, the current I k carried by the electrode E k can also be split 
into two parts, I k {t) = i k {t) + 1% (t). The first term is referred to as the inner current, 
the second is the vacuum part. The vacuum part is nothing but the displacement current 
related to the vacuum field <|>( vac ). It is shown in the appendix that it has a purely capacitive 
character. One obtains 

1=0 

The coefficients C k \ are the capacitance coefficients defined by 



C kl = / # fc e e r WrdV. (20) 

JdV 

It is also be shown in the appendix that the vacuum current satisfies Kirchhoff 's law of 
current conservation, i.e., the sum of the currents over all electrodes (including grounded 
electrodes) is equal to zero, 



With the help of the characteristic functions introduced in the previous section, the inner 
part of the current ij~ can be written as a flux integral over the electrode 



i k = - V^ k -jd 3 r. (22) 



V 

(The explicit derivation can be found in the appendix.) The inner part of the current obeys 
Kirchhoff's law of current conservation, i.e., their sum over all electrodes (including ground 
electrodes) is equal to zero (see also appendix). 

We now establish energy balance equations for the system. Owing to the fact that the 
coupling of the vacuum field and the inner potential vanishes (see appendix), we can split 
the electrical energy of the system into a vacuum part £'( vac ) and an inner part given by 

e ( el ) = ! I eo£r V0-V0rf 3 r (23) 
Jv 2 

and 

1 N N 

Eivac) = _j2Y^C kl V k V l . (24) 

k=0 1=0 

The vacuum part can be expressed in terms of the previously defined capacitance coefficients 
(12"U]I . For the vacuum part, the energy balance simply describes the charging and decharging 
of the linear multiport capacitor Cm, 

jp(vac) N N RV N 

k=0 1=0 1=0 

The inner energy e is the sum of the inner electrical energy and the kinetic energy due 
to the motion of electrons 



£ = + e (kin) = / eoerV 0. V rf 3 r+ / -J 2 d 3 r. (26) 

Jv 2 Jp 2e o^ c 

The temporal change in the inner energy is governed by the balance of the power input and 
the Ohmic dissipation due to collisions 

S = X>- / -^rfdh. (27) 
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IV. A LUMPED ELEMENT ANALOGY 



As a simple problem and to introduce an abstract mathematical formulation, we consider 
the response of a lumped series resonance circuit driven by a voltage V(t) as shown in figure 
HJ In terms of active resonance spectroscopy one can interpret this scenario as a one-port 
concept which does not couple to the outer ground. We assume that only one voltage source 
and one resonance mode is active. For the above mentioned exemplary character we assume 
the values of the lumped circuit elements as given. (To evaluate real plasma diagnostic 
concepts the values have of course to be determined by means of a mathematical model.) 

The exemplary model can be cast as a system of two coupled ordinary differential equa- 
tions of first order for the capacitor voltage U and the current /, 

H-c> (28) 

dI U RI l t ,M 

Tt = "I " T + I v{t) - (29 » 

To solve the system the Fourier transform can be used and one obtains its admittance that 
is interpretable as the system response of the one-port probe 

Y{u) = ^ = (30) 

R + iuL + — — 

The same result is deriveable by an abstract mathematical formulation of the system. To 
introduce this formulation we refer to the ordered pair of U and I as the state of the system, 
denoted by the state vector 

\z) = \u,i)= r j . (3i) 

Similarly, we introduce the vector of the excitation of the system 



e) = i . (32) 




The operators of the non-dissipative and the dissipative dynamics are given by 



c 





Td = ( q _R | - (33) 
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With the help of these definitions, the behavior of the lumped element circuit can be ex- 
pressed in terms of a compact dynamic equation, 

^=T c \z}+T D \z) + V(t)\e}. (34) 

Since the set of all states \z) forms a vector space Z it is advantageous to employ as the 
underlying scalar field the set of complex numbers C. To turn Z into a Hilbert space, the 
appropriate scalar product has to be defined. Due to the fact that the dimension of the 
vectors \z) is inhomogeneous, the Euclidian product is not available. However, a physically 
reasonable quadratic form exists, i.e., the energy e = (CU 2 + LI 2 ) /2, which is stored in the 
capacitance and the inductance of the resonant circuit. This motivates the definition of a 
scalar product of two state vectors \z) and \z) as 

(z\z) = (U, I\U, I) := CU*U + LI* I. (35) 

The energy is a non-negative real quantity even for complex states, 

e = \M\ 2 = \{z\z). (36) 

The defined scalar product is compatible with the dynamical equations. It can be shown 
that the operator Tc is anti-Hermitian with respect to the scalar product, 

(T c z\z) = -(~z\T c z), (37) 

while the dissipative operator To is Hermitian and negative definite, 

(T D z\z) = (z\T D \z), (38) 
(z\T c \z) < 0. (39) 

Furthermore, it can be shown that the observable response of the system, i.e., the current /, 
is the scalar product of the excitation vector and the state vector \z). Thus, the observation 
vector is equal to the excitation vector, 

I=(0,^\UJ) = (e\z). (40) 

The energy balance equation ( 13"B"|) can be transformed to the well known energy balance 
equation of a series resonance circuit 

Jt( 1 2 CU2+1 2 Ll2 )=- Rl2 + IV{t) - (41) 



To analyze the response of the system we apply the Fourier transform to (134]) . The 
current is then given by 

I = (e\z) = (eh ^—r\e)V(uj) = Y(u)V. (42) 

iuj - T c - T D 

The response function Y(u) of the system can be directly specified by evaluating the matrix 
products. It should be noted, that two extra matrices from the scalar product (1351) appear 
in the calculation. However, the response function can also be calculated by expanding the 
state vectors into the eigenfunctions of the operator. One ultimately obtains exactly the 
same expresion than in (130]) . This result of the simple RCL-series circuit allows the physical 
interpretation of the abstract solution that is derived in a section further down. 



V. FUNCTIONAL ANALYTIC FORMULATION OF THE PROBE RESPONSE 

We are now equipped with all the material necessary for a functional analytic formulation 
of the system of dynamical equations (TTTT) and (Ti~8l) . We interpret the charge density p(r) 
and the current density j(f) defined on the plasma domain P as the variables of the system 
state \z). (It is advantageous to allow for complex states; of course, physically observable 
states are always real.) In Dirac's bra-ket notation we define an appropriate function space 
as 

H={\z) = \pJ) = (p,j) T }. (43) 

Using two distinct system states \z) = \p,j) and \z') = \p',j') we define a scalar product 
{z\z') which is motivated by the inner energy functional ( 126]) . The following Hermitian 
form is compatible with the dynamical equations and satisfies all requirements of an inner 
product, namely i) conjugate symmetry, ii) sesquilinearity, and iii) positive definiteness. 

(z\z')= [ eoer V0{p}*-V0{p / }rf 3 r+ / -l^f.j'd 3 r. (44) 
Jv Jp eo^pc 

The scalar product defines a corresponding norm ||z||. The norm-square of a state is equal 
to its inner energy, up to a factor of two 

\\z\\ 2 = (z\z) = [ eoer V0{p}*-V0{p}rf 3 r+ / -^f.jd 3 r = 2e. (45) 

JV JP e 0^ pc 

Equipped with this scalar product, and suitably completed, the defined function space of 
all possible states % satisfies the properties of a Hilbert space H. It should be noted that 
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the elements of % are not all regular functions. The charge part may contain also singular 
functions like surface charges. However, this acceptable as long as the norm ||z|| is finite. 

One important family of states \e k ) can be obtained from the characteristic vacuum 
functions ^ k which have only a current part, their charge part and also the corresponding 
potential are equal to zero, 

|e fc ) = |0,-e ^ e W fc ). (46) 

The scalar product of je^) with a system state \z) is an expression for the inner current i k 
which is carried by the electrode E k . Utilizing the definition of the scalar product fj4"5]) one 
obtains 

(e k \z) = [ W k -jd 3 r = i k . (47) 



The dynamical equations ( TT71) and ( 1T81) can be written in a compact form [23[. Therefore, 
we define two dynamical operators Tc and Tp which act on the Hilbert space H. The 
conservative operator Tc contains the effects of the electron inertia and the electric field; it 
is defined by the action 

Tc\p,l) = I - V-J, -eoa&WM). (48) 

Tc can be shown to be anti-Hermitian, i.e., Tc = —T c (see appendix). Thus the associated 
eigenvalues are purely imaginary and represent the frequency behavior. The dissipative 
operator Tjj contains the effects of the collisions 

T D \pJ) = 10,-i/J). (49) 

It can be shown that Td is Hermitian, i.e., Tjj = T^ D (see appendix). The eigenvalues are 
negative definite, (z\T^\z) < 0. Furthermore, they are real and contain information about 
damping of the system. With the help of these definitions, the dynamical equations (JTJ 
and ( 1181) of the system can be written as 

= J-C\Z] T J-D\Z] -T 

1=0 

Applying a Fourier transform to this abstract form and solving for \z), the current as the 
desired measure can be obtained 



T c \z)+T D \z) + J2Vih)- (50) 



N N 

ik = (e k \z) 



i=i 



ii 



Thus, the response function, can be written as 

Um(uj) = (e k \ (iu - T c - To)' 1 \ei). (52) 

The main result of this analysis is that the response function yki{oo) is given by the 
matrix elements of the resolvent of the dynamical operator equation evaluated for values 
on the imaginary axis. The result can be ultimately interpreted using equivalent lumped 
element circuits as described above. 



12 



VI. RESULTS FOR A ONE-PORT PROBE 



With this result and the insights gained from the simple application of the formulated 
theory, we are able to discuss the result of the abstract model, given by equation (J5TJ) in 
terms of a lumped element circuit model. For a one-port probe (we account for two electrodes 
plus ground) the current to one of the two powered electrodes is given as 

2 2 

ik = {e k \z) = ) y {e k \ (ico - T c -Toy 1 \ei)Vi = ^ y k iVi- (53) 

1=0 1=0 

Although the response function yki{u) has the same structure than (|42l) . It represents the 
coupling between the powered electrodes directly and the coupling from each electrode to 
ground. In this case, the dynamical operator T = Tq + To is not a matrix as in the simple 
example, but a functional operator. The current is able to flow on infinite different paths 
through the plasma, each path can be expressed by a series resonance circuit. Therefore, 
the response function represents an infinite number of series resonance circuits connected in 
parallel. 

In summary, the complete equivalent circuit consists of a branch which represents the 
direct coupling between the two powered electrodes. Two additional branches mimic the 
parasitic coupling of the electrodes to ground or "infinity" . Each branch itself is complex, 
and the relation between the plasma and probe parameters and the values of the circuit 
elements are generally very complicated. This case is depcited in figure [3j For realizations 
with more than one port the coupling gets even more complex. 

Under strictly electrical and geometrical symmetry the situation for a one-port probe 
becomes simpler. It can be shown via a A-Y transform (circuit theory) that the coupling 
to ground or "infinity" vanishes. As the field decays quickly with increasing distance from 
the probe, the measurements become actually local. The equivalent circuit becomes more 
simple: an infinite number of series resonance circuits, each of which representing a discrete 
resonance mode, are connected in parallel to a vacuum coupling. 

The multipole resonance probe proposed by the authors 17| fulfills these requirements to 
a great extent - it consists of two separated hemispherical electrodes covered by a dielectric 
shielding. Thus its geometry is approximately symmetric (with respect to rotations around 
the center of the probe) and its electrical behavior is completely symmetric (with respect to 
the mapping ujt — > ut + n and $ — > — $. Additionally, the spherical geometry allows for a 
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compact mathematical description of the resonance behavior and simple algebraic expression 
for its resonance frequencies. 



VII. SUMMARY / CONCLUSION / OUTLOOK 

We have formulated a model of electrostatic active resonance spectroscopic methods in 
terms of a functional analytic approach. The abstract character of these methods allow for 
a generalized model of the different realizations as the geometrical and electrical character 
do not enter explicitly in the model equations. Instead of a numerical solution we were able 
to present a general analytic solution for the system state that gives physical insight in the 
resonance behavior. The scalar product between this state and the observation state was 
shown to be the inner current at one electrode. Part of the current is the inner admittance 
that is determined as a matrix element by the resolvent and the observation states. This 
result can be interpreted in terms of a lumped element equivalent circuit. To guide the 
physical interpretation of the functional analytic approach this method was explicitly used 
to obtain the systems response of a idealized situation, i.e., that the equivalent circuit is 
described by a simple RCL-series circuit. 

As an example we presented the equivalent circuit of an one-port probe. An analysis of 
this circuit allows for an optimization of the probe geometry. In terms of such an optimiza- 
tion, a coupling from the electrodes to the outer ground is obstructive, hence a optimized 
probe design has to prevent this coupling. Due to a symmetrical geometry and a symmetrical 
applied voltage it is possible to transform the circuit in an equivalent circuit of one branch 
between the two electrodes. This branch is given by an infinite number of parallel series 
resonance circuits and each represents a resonance mode parallel to a vacuum coupling. The 
authors proposed the multipole resonance probe [l^l as an optimized concept. 

Measurements with the prototype of the multipole resonance probe 16[ have shown the 
feasibility of this concept. In further measurements to determine the collision frequency of 
electron-neutral collisions in low pressures we observed a stronger damping of the resonance, 
which means a higher collision frequency, than the model predicts. This demonstrates 
that kinetic effects play a crucial role on the damping of the resonances in low preassures. 
Here, only a kinetic treatment of the electron dynamics can describe the resonance behavior 
correctly. First attempts to establish a kinetic theory based on functional analytic methods 
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are in progress. 



Appendix A: Derivation of the balance equations 

In the following, the balance equations and their relations will be formally derived. 
The calculations utilize the facts that the particle current vanishes at the electrodes 
due to the dielectric shielding, the total inner current is divergence-free, the fy k are the 
characteristic functions defined by eq. ( fTBl and employs the boundary conditions of $ and <p. 



The vacuum current can be expressed via capacitance coefficients in the following form 

,r> . -/ -i*±% -f>f . (A1) 

It can also be shown, that Kirchhoff 's law of current conservation is conserved. 

TV N N 



J2°ki = I y k eoe r Wrd 2 r = J2 [ V-(* fc e e r W,y 

k=0 k=0 JdV k=0 Jv 

N . N . 

= Y] / e e r W fc -V* z d 3 r = V / e e r W fe * r <?', 

u-n JV u _ n JdV 



k=0 " y k=0 



e e r V^ r d 2 f= / V- (e e r Vty) d 6 r = 0. (A2) 
av Jv 

The concise for the inner current can be derived as following 

f <9V</> f (-; 0V<A i* (-t dV<f)\ l3 

.* = - eoer-Qj-d r = -J dv * fc - ^— J -d r = - ^ V**- [j - e e r — j d r 

^ k -jd 3 r + ^W fc -Ue r ^j d 2 f=- J W k -jd 3 r. (A3) 

Also the inner current obeys Kirchhoff's law 

It can be shown, that the coupling of vaccum- and inner energy vanishes 

E (cpl) = j e o e r V0-V$ (vac) d 3 r 
Jv 

= [ 0e o e r V$ (vac) -d 2 f- I </>V-(e e r V$ (vac) )d 3 r = 0. 

JdV Jv 
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(A5) 



For the inner energy, we obtain after some calculation the following relation which reflects 
the Ohmic dissipation due to the collisions: 



e o e r V0-— -^tfr + / Y3'it dr 

at Jv at Jp e u£ c at 



Jp V / 

W5 + '^ + j[( j -^){-H A -/ 



TV 



AT 

V 



= E^- / -V?^- (A6) 
Appendix B: Properties of the dynamical operators 

In this section the properties of the dynamical operators are derived. It can be shown, 
that the operator of the conservative dynamics is anti-Hermitian 



(z'\T c \z} + (z'\Tl\z}= ! eoer V0{p'}*-V0{-V-jVr+ / -L-j'* ■ (-e o c^ e V0{p})d 3 r 

Jv JP € 0U pe 

+ / eo e r V0{-V-j / }*-V0{p}d 3 r+ f _i^(- eo ^ e V0{p'})*-Jd 3 r 
Jv Jp e oU pc 

= f rtp'YV- (eoe r V0{V-J}) tfV + Jv-j'^p^r 

+ Jv- (eoe r V0{V-J'}*) 0{p}d 3 r + J<f){p'yV -jd 3 r = 0, (Bl) 
and the operator of dissipative dynamics is Hermitian 

(z'\T D \z) - (z'\ T% \z) = [ -^j'*-(-vj)d 3 r - f -L- (-uj'*) -j d 3 r = (B2) 

JP^pe V J JP^pe V 7 

Finally, we can show that it is negative definite 

(z\T D \z) = [ -\f-(-uj)d 3 r < 0. (B3) 

Jp e W pe V 7 
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FIG. 1: Illustration of the abstract model for a N-electrode system. The electrodes are shielded to 
each other and the plasma by some dielectric medium. The whole probe, surrounded by a plasma 
sheath is immersed in a plasma volume. 
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FIG. 2: A lumped series resonance circuit driven by a voltage V{t). This circuit is interpretable as 
an one-port concept of active resonance spectroscopy which does not couple to the outer ground. 
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FIG. 3: Equivalent circuit for an one-port probe. It consists of a branch which represents the direct 
coupling between the two electrodes, and two additional branches which stand for the parasitic 
coupling of the electrodes to ground or "infinity". Each branch in itself may be complex, and 
the relation between the plasma and probe parameters and the values of the circuit elements is 
generally very complicated. 
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